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Abstract 

Within the context of linear elasticity we show that a two-terminal network of 
springs and masses, can respond exactly the same as a normal spring, but with a 
frequency dependent spring constant. As a consequence a network of such springs 
can have a frequency dependent effective elasticity tensor but zero effective mass 
density. The internal masses influence the elasticity tensor, but do not contribute 
to the effective mass density at any frequency. 

Keywords: metamaterials, effective mass density. 

Introduction 

The mass density term which enters the macroscopic wave equation of linear continuum 
elastodynamics need not be the same as the mass density calculated by local volume 
averaging the interior masses. This was observed by Berryman [1] in the context of an 
approximation scheme for suspensions of particles in an inviscid fluid and later found to be 
more generally true for composites constructed from components having a high contrast 
in stiffness and density. The reason is simply that interior masses do not necessarily 
locally move together in lock step motion with each other, even when the wavelength 
is long, and therefore the locally averaged momentum is not simply a product of the 
locally averaged mass times the locally averaged velocity [2j. At a given frequency, the 
effective mass density entering the wave equation can be anisotropic, negative, or even 
complex [31 SJ [5j El [?1 [2i [8] . I n general, the vibrations of interior masses cause the effective 
elasticity tensor and the effective mass density to depend on frequency. This raises the 
interesting question, which will be explored here, as to whether internal masses can 
influence the effective elasticity tensor, and cause it to be frequency dependent, but not 
contribute to the overall effective mass density at any frequency? We will see the answer 
for mass-spring networks is yes, theoretically it can, provided we are strictly working 
within the framework of linear elasticity and ignore gravity and non-linear instabilities 
such as buckling. 
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Dispersive Normal Springs 



If we have a normal spring with spring constant k and terminals at the points xi and X2 
then the response of the spring, relating the two forces fi and f 2 applied to the spring at 
the terminals to the displacements Ui and u 2 there, takes the form 



kn ■ (ui — u 2 ) 
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where n 



xi - x 2 
|xi - x 2 | 



(1) 



We will define a dispersive normal spring to be any two terminal network of springs and 
masses such that the relation between the forces and displacements at the terminals at 
any frequency takes the form (JT]), but with k being dependent on frequency. In particular, 
the relation ([!]) implies that the forces fi and f 2 are equal in magnitude and opposite in 
direction, and directed parallel to xi — x 2 . 

An essentially explicit scheme for constructing dispersive normal springs is a corollary 
of Theorem 4 of Guevara Vasquez, Milton and Onofrei [9J. That theorem covers the much 
more general case of characterizing all possible time-dependent responses of multiterminal 
mass-spring networks, and generalized earlier work of Camar-Eddine and Seppecher [10] 
on characterizing the response of multiterminal spring networks. Our objective here is 
to present a simpler construction where the mechanism responsible for the behavior of 
the network is easy to grasp. 

To construct a dispersive normal spring, consider first the four terminal spring network 
of figure 1, consisting of a square diamond of identical springs with four other identical 
springs as legs, which extend directly outwards from the four vertices of the diamond to 
the terminals. This network is a rank one network: it only supports one loading (and 
all multiples of this loading). To see this suppose the spring linking terminal 1 to the 
diamond is under tension T. Then, by balance of forces, the four springs in the diamond 
must be under tension T / v^2 and the other legs joining the diamond to the terminal edges 
must be under tension T. Thus forces at all the four terminal nodes are determined if 
we know the force at one terminal node : the response of this network, relating the 
four forces fi, f 2 , f3 and f 4 at the terminals to the displacements Ui, u 2 , u 3 and U4 
there necessarily takes the form (f 1; f 2 , f 3 , f 4 ) = (Tn, — Tn, Tn^, — Tn±) where the scalar 
tension T depends linearly on ui, u 2 , 113 and u 4 (here n denotes the unit vector pointing 
from x 2 to xi and ni is the orthogonal unit vector pointing from X4 to X3). Thus the 
response matrix W which relates the forces at the terminals to the displacements there, 
via the linear relation, 



(2) 



is a rank-one matrix. Since (no matter what the network), this response matrix is a 
positive semidefinite symmetric matrix, we deduce that 
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where the constant g > scales in proportion to the stiffness of the springs in the 
network. Now, let us place a mass m at the terminals 3 and 4, and make them interior 
nodes, and consider the response at frequency u of the resulting two terminal network. 
The only forces at nodes 3 and 4 are the inertial forces f3 = mw 2 U3 and f*4 = mw 2 U4. 
Substituting this in (j3]) and eliminating u 3 and u 4 from the resulting equations, gives a 
relation between (f l5 f 2 ) and (ui, u 2 ) exactly of the form (JTJ with a spring stiffness 

2gmu 2 

k{u) = — — . 4 

mu z — 2g 

Thus this dispersive normal spring behaves exactly like a normal spring but with a 
frequency dependent spring stiffness, which will be negative for frequencies below uq = 
■sj2gjm. The masses in the dispersive normal spring do not move, to first order in the 
displacement, when the two terminal nodes undergo a rigid body motion. They only 
move when the spring is extended or compressed. 

The dispersive normal spring constructed here is floppy. For example, within the 
framework of linear elasticity, the interior diamond can be infinitesimally rotated with 
no change in the forces fi, f 2 , f3 and f 4 and displacements Ui, u 2 , 113 and 114 at the 
terminals. Strictly speaking one should perturb the network by adding a scaffolding of 
additional springs with very small spring constants to uniquely determine the interior 
displacements, but we refrain from doing so as to avoid complications. (Since then fi will 
not be exactly equal to — f 2 .) The undetermined displacements do not effect the overall 
response of the dispersive normal spring. However associated with this degree of freedom 
is a buckling mode: when the springs in the network are under compression rather than 
tension one can expect, within the framework of non-linear elasticity, that the interior 
diamond will rotate, one way or the other, to relieve this compression. We ignore this 
since we are working only within the framework of linear elasticity. 

We will now forget about the internal structure of the dispersive normal spring and 
treat it as a single object. 



A material having frequency dependent elasticity ten- 
sors and zero effective density 

Suppose we have have a periodic network (such as a triangular network) of identical 
normal springs with spring stiffness k. Let the nodes of this network be called "primary 
nodes" . The effective elasticity tensor of the network C will be proportional to k and we 
can write C = kC°. Now if we replace each spring in the network by a dispersive normal 
spring with the same spring stiffness k(cu) given by (J4]), then the resulting material will 
have a frequency dependent effective elasticity tensor 

C = k{u)C° = 2gmu 2 C°/{mu 2 - 2g), (5) 

which is negative definite for frequencies below u = ^J2g/m. 

At the same time, within the framework of linear elasticity, the effective density of 
this network will be zero! By this we mean that at any fixed frequency, and to first order 
in the displacements, that the macroscopic response of the network will be the same as 
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Figure 1: Sketch of the dispersive normal spring. The open circles represent terminal 
nodes, and the closed circles could be either terminal nodes or interior nodes with masses 
attached. The straight lines represent springs. The large arrows represent external or 
inertial forces acting on the nodes at one instant in time. The two small arrows on each 
spring give the direction of the force which the spring exerts on the node nearest to the 
arrow. 



that of a periodic network of identical normal springs, without masses at the nodes. The 
internal masses cause the effective elasticity tensor to depend on frequency (and cause it 
to be negative definite for u < u> ) but do not contribute to the effective density, since 
the internal masses do not move (to first order in the displacements) when the primary 
node lattice is translated. Of course, as opposed to normal linear elasticity which can 
be useful even if there are large displacements such as rotations, the linear elasticity 
approximation here will only be valid for displacements which are small compared to the 
length of the springs in the network, if at all. 

Some care needs to be applied in this notion of effective elasticity tensors. If we take a 
large sample of the periodic network of dispersive normal springs then we need to ensure 
that any cut dispersive normal spring at the boundary of the sample is removed and that 
there are no interior nodes of the dispersive normal springs in contact with the surface 
loadings applied to the boundary of the sample. If there are body forces (perhaps due to 
electric field gradients and polarizable nodes) then we need to make sure these only act 
on the primary nodes common to the original network of normal springs and not on the 
internal nodes and masses of each individual dispersive normal spring. 
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